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Abstract 

Transport equations can be derived from quantum field theory 
assuming a loss of information about the details of the initial state 
and a gradient expansion. While the latter can be systematically 
improved, the assumption about a memory loss is in general not 
controlled by a small expansion parameter. We determine the range 
of validity of transport equations for the example of a scalar 
theory. We solve the nonequilibrium time evolution using the three- 
loop 2PI effective action. The approximation includes off-shell and 
memory effects and assumes no gradient expansion. This is compared 
to transport equations to lowest order (LO) and beyond (NLO). We 
find that the earliest time for the validity of transport equations is set 
by the characteristic relaxation time scale tdamp = — where 

—T!f^/2 denotes the on-shell imaginary-part of the self-energy. This 
time scale agrees with the characteristic time for partial memory loss, 
but is much shorter than thermal equilibration times. For times larger 
than about tdamp the gradient expansion to NLO is found to describe 
the “full” results rather well for < 1. 
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1 Introduction 


Important phenomena in high-energy physics related to collision experiments 
of heavy nuclei, early universe cosmology and other complex many-body 
systems urge a quantitative understanding of nonequilibrium dynamics in 
quantum held theories. One of the most crucial aspects concerns the 
characteristic time scales on which thermal equilibrium is approached. Much 
of the recent interest derives from observations in collision experiments of 
heavy nuclei at RHIC. The experiments seem to indicate early thermalization, 
whereas the present theoretical understanding of QCD suggests a signihcantly 
longer thermal equilibration time pQ. 

Most theoretical estimates concerning the important question of 
thermalization have been obtained from transport theory, or assuming in 
addition a quasi-particle picture, from kinetic theory. Transport equations 
can be derived from quantum held theory assuming a loss of information 
about the details of the initial state and a gradient expansion. While the 
latter can be systematically improved, the assumption about a memory loss 
is in general not controlled by a small expansion parameter. The time for 
a (partial) loss of memory represents a minimum characteristic time after 
which transport descriptions may be applied. It is important to quantify 
this characteristic time in quantum held theory. The question is whether 
transport or kinetic theory can be used to quantitatively describe the early- 
time behavior, which is necessary if their application to the problem of fast 
apparent thermalization is viable. 

With the advent of new computational techniques a direct account of 
quantum held degrees of freedom becomes increasingly feasible. There has 
been important progress in our understanding of nonequilibrium quantum 
helds using suitable resummation techniques based on 2PI generating 
functionals [2]. They have been successfully used to describe far-from- 
equilibrium dynamics and subsequent thermalization in quantum held 
theory muiainiiaiHi. While it is difficult to test analytic approximations 
for theories such as QCD, where strong interactions play an important role, 
the description of scalar and fermionic theories for not too strong couplings 
seems well under quantitative control by now. 

In this work we determine the range of applicability of transport equations 
for a scalar 5 '^$^ theory. The latter has been frequently studied in the 
past using the required 2PI loop- jU] or 2PI 1/^expansions beyond lowest 
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order laiini.' Here we employ a three-loop 2PI effective action in 3 -|- 1 
dimensions. It includes direct scatterings, off-shell and memory effects. 
Most importantly in this context, it employs no derivative expansion or 
quasi-particle assumption. We solve the nonequilibrium evolution for a 
class of nonequilibrium initial conditions reminiscent of some aspects of the 
anisotropic initial stage in the central region of two colliding wave packets. 

We then derive the corresponding transport equations to lowest order 
(LO) and next-to-lowest order (NLO) following standard prescriptions. Both 
LO and NLO are then compared to the “full” results as described in Sec. HI 
It turns out that the NLO corrections are crucial already for relatively small 
couplings, which made it necessary to go beyond the LO study presented in 
Ref. m The LO expressions are typically employed to obtain kinetic or 
Boltzmann equations using in addition a quasi-particle approximation. We 
emphasize that our comparisons do not rely on any kind of quasi-particle 
ansatz. We concentrate on the applicability of transport theory based on 
a derivative expansion, which also provides the necessary condition for the 
applicability of kinetic equations. 

2 Exact evolution equations 

We start by considering the exact nonequilibrium evolution equation for the 
time-ordered two-point correlation function 

G{x,y) = (T$(x)$(i/)) 

= F{x, y) --p{x,y) sign{x°-y^), (2.1) 

i.e. the expectation value of the product of two Heisenberg field operators <I> 
with time-ordering T. Here (...) includes the trace over the density matrix 
describing the initial state j2]. For simplicity, we will consider a scalar 
field theory in the symmetric regime where the field expectation value 
vanishes, i.e. (‘F(x)) = 0. 

In general, there are two independent two-point functions, which can 
be associated to the commutator and anti-commutator of two fields. 
Accordingly, the second line in Eq. (EH) is a rewriting, which expresses the 
time-ordered product using 

_ P{x,y) = i([$(x), $(?/)]), (2.2) 

^Here N denotes the number of field components. 
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F{x,y) 




(2.3) 


Here p(a;, y) denotes the spectral function and F{x, y) the statistical two- 
point function. For the real scalar theory these determine the imaginary 
part and the real part of the two-point function (EH), respectively. While 
the spectral function encodes the spectrum of the theory, the statistical 
propagator gives information about occupation numbers. Loosely speaking, 
the decomposition EH makes explicit what states are available and how 
often they are occupied. For nonequilibrium, F{x,y) and p{x,y) are in 
general two independent real-valued two-point functions. A simplihcation 
occurs in vacuum or thermal equilibrium where F{x, y) and p{x, y) are related 
by the fluctuation-dissipation relation [B]. We note that the spectral function 
encodes the equal-time commutation relations: 

2 /) Uo=yO = 0 , d^op{x,y%o=yO = 5{x-y). (2.4) 

The exact evolution equations for the statistical and spectral functions 
read (for a detailed derivation see e.g. [2]): 

poo 

[D^ +M‘^{x)] F{x,y) = - dz° / d'^zJ:p{x,z)F{z,y) 

Jo J —oo 

p'jP poo 

+ dz^ d‘^zJ:p{x,z)p{z,y), (2.5) 

Jo J —OO 

px^ poo 

[D^ +M‘^{x)] p{x,y) = - dz^ d‘^zEp{x, z)p{z,y). ( 2 . 6 ) 

J J—oo 

These are causal equations with characteristic “memory” integrals, which 
integrate over the time history of the evolution starting at time to = 0 . 
Since they are exact they are equivalent to any kind of identity for the two- 
point functions such as Schwinger-Dyson/Kadanoff-Baym equations. Eqs. 
(EH and (EH are obtained for Gaussian initial conditions, which are also 
underlying transport equations (cf. Sec. Ej)- The restriction to a class of 
initial conditions itself represents no approximation for the dynamics and 
for non-zero self-energies higher irreducible correlations build up for times 
t > 0.2 

^More involved initial conditions can be considered using higher n-PI effective actions, 
which would involve evolution equations for higher n-point functions beyond the two-point 
correlator laiiaiinj. 
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The spectral (Sp) and statistical self-energy parts (Sf) are obtained 
from the proper self-energy E(x, ?/), which s um s all one-particle irredncible 
diagrams, using a decomposition similar to Eq. (EH) 0: 

T,{x,y) = —iTS^\x)5{x — y) 

+T.f{x, y) - ^Ep(x, y) sign(x° - ?/°). (2.7) 

Since corresponds to a space-time dependent mass-shift it is convenient 
to introduce the notation 

M^{x]G)=m^+ T.^^\x]G). ( 2 . 8 ) 

To close the set of equations of motion (EH and (EH we need further 
equations that connect the self-energies to F and p. This link is established 
by defining an approximation scheme based on the two-particle irreducible 
(2PI) effective action jH],^ which to three-loop order in the scalar theory 
yields the self-energies 

Sf(x,2 /) = -QQg^F{s,x) (^‘^{x,y)-^p^{x,y'^ , (2.9) 

^p{x,y) = -2SSg^p{x,y)(^‘^{x,y) - ^p^{x,y'^ . ( 2 . 10 ) 

This approximation has been shown to describe thermalization ininiEi and 
includes direct scatterings, off-shell and memory effects. Therefore, it goes 
far beyond standard kinetic descriptions or Boltzmann equations. Most 
importantly in this context, it employs no derivative expansion. The latter 
is a basic ingredient for transport or kinetic theory. 

^For solving Eqs. EH-EH we introduce a further approximation by restricting the 
memory time integration to a maximum interval in order to fit in the memory of the 
computing device. In this paper we use a 8192^ x 16^ lattice using about 100 GB on 
a distributed memory super-computer. The used integration interval is ~ 1200 to^^ for 
= 0.5 and ~ 1700m^^ for = 1. This interval is much longer than the characteristic 
memory or damping time of the propagator, thus this memory-cut has a negligible effect 
on the dynamics. 
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3 Transport equations 

The spectral function (Q is directly related to the retarded propagator, 
Gr, or the advanced one, Ga, by 

GR{x,y) = -y^)p{x,y), GA{,x,y) = -Q{y^-x^)p{x,y). (3.1) 

Similarly, the retarded and advanced self-energies are 

^R{x,y) = Q{x^-y^)T.p{x,y), T.A{x,y) =-Q{y^ - x^)T.p{x,y). (3.2) 

With the help of this notation, interchanging x and y in the evolution 
equation (ESI) and subtraction one obtains 

(□, -Uy + M^ {x) - m 2 {y) ) F {x, y) 

= j d'^+^ 2 ; 0 ( 2 ;°) [f {x, z) Sa {z, y) -h Gr {x, z) 'Lr {z, y) 

- T,r {x, z) F (z, y) - T,f {x, z) Ga {z, y) j . (3.3) 

The same procedure yields for the spectral function 

- Dy + m 2 (a:) - m 2 {y)^p{x, y) 

= J d ‘^+^ 2r (Gr {x, z) Ep {z, y)+ p {x, z) Ea {z, y) 

- Ep (x, z) Ga {z, y) - Er (x, z) p {z, y) j . (3.4) 

So far, the equations (ESD and (El are fully equivalent to the exact equations 

(El and (ESI). 

3.1 Assumptions 

Transport equations are obtained from the exact equations (El and (El 
by the following prescription: 

1) Drop the 6'-function in Eq. (13.3j) . which sets the time to = 0 where the 
system is initialized. Dropping it amounts to sending the initial time 
to the inhnite past, i.e. to —> — 00 . Of course, a system that thermalizes 
would have reached already equilibrium at any hnite time if initialized 
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in the remote past. Therefore, in practice a “hybrid” description is 
employed: transport equations are initialized by prescribing F, p and 
derivatives at a finite time using the equations with to as an 

approximate description. 


2) Employ a gradient expansion to Eqs. (Id.d|l and (ld.4jl . In practice, this 
expansion is carried out to lowest order (LO), i.e. no derivatives, or 
next-to-lowest order (NLO) in the number of derivatives with respect 
to the center coordinates 


= 


2 


and powers of the relative coordinates 


(3.5) 


st‘ = x“-yr 


(3.6) 


3) Even for hnite X° one assumes that the relative-time coordinate 
ranges from —oo to oo in order to achieve a convenient description 
in Wigner space, i.e. in Fourier space with respect to the relative 
coordinates (jSSD. 

We emphasize that the ad hoc approximations 1) and 3) are in general 
not controlled by a small expansion parameter. They require a loss of 
information about the details of the initial state. More precisely, they can 
only be expected to be valid for sufficiently late times times t when initial¬ 
time correlations become negligible, i.e. (<h(0,x)<h(f,y)) ~ 0. The standard 
approximations 1) and 3) may, in principle, be evaded. However, if they are 
not applied then a gradient expansion would become too cumbersome to be 
of use in practical calculations. 

The derivation of transport equations has been discussed in great detail 

in the literature [HI UHl 1^ El El El 120112111221120112111201 120] • Most 

discussions focus on kinetic theory employing the additional approximation 
of a suitable quasi-particle ansatz, which goes beyond assumptions 1) - 3). 
In order to be general, we do not restrict the following analysis to any quasi¬ 
particle picture. 

For the description in Wigner space we introduce the Fourier transforms 
with respect to the relative coordinates, such as 

/ 2X° 1-00 

ds°e*“^"° / d'^s e*P*F(X +s/2, X-s/2) 

-2X0 J_ao 
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Figure 1: For given finite 2X^ = + ip the relative coordinate has 

a finite range from —2X^ to 2X^. 


f2X° 


= 2 


ds° cos(o^s°) / d'^s PP^F{X + s/2, X - s/2), (3.7) 


using the symmetry property F{x,y) = F{y,x) for the second line. We 
emphasize that the time integral over s° is bounded by The time 

evolution equations are initialized at time x^ = y^ = 0 such that > 0 and 
> 0. According to Eq. (uni), the minimum value of the relative coordinate 
is then given by —y^ = —2X° for = 0 while its maximum value is 
= 2X‘’ for = 0 as illustrated in Fig. ^ Similarly, we dehne 


/ 2X^ poo 

ds^F^^° / d‘^sFP^p{X + s/2,X-s/2) 

-2X0 
1-2X0 

Jo 


ds° sin(a;s°) / d'^s+ s/2,X — s/2), (3.8) 


where a factor of i is included in the dehnition to have g real and we have 
used p{x,y) = —p{y,x). It will be convenient to also introduce the cosine 
transform 

p2X^ poo 

g+{X;uj,p} = 2 / ds° cos(a;s°) / d'^se*^^p(X + s/2,X — s/2). (3.9) 

Jo J—oo 


The equivalent transformations are done to obtain the self-energies 
Ef-(X;u;,p), E^(X;uJ,p) and the cosine transform E^(X;uJ,p) corresponding 
to Eq. (13.9j) . From Eqs. dSH) and dH it follows that the cosine transforms 
are related to the real part of the retarded functions 


p+(X;a;,^ = 2ReGR(X;a;,^ , S+(X;a;,^ = 2ReSi^(X;a;,^ . (3.10) 
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In order to exploit the convenient properties of a Fonrier transform it is 
a standard procedure to extend the limits of the integrals over the relative 
time coordinate in Eqs. (HID, (EHl) and (jH to ±cxo. For instance, using the 
chain rule 

□ 3 ; - Dy = 2ds^dxi^ (3.11) 

and the gradient expansion of the mass terms in Eq. (Id.d|l to^ NLO: 

(X - 0 ~ (X) , (3.12) 

then the LHS of Eq. (13.3|) becomes in Wigner space 


drf+ig ^ips ^2ds^dx^ + s>^dx>^M^ (X)] F (X + s/2, X - s/2) 

) 

= -i + {dxt^M‘^ (X)) dp^] F{X-,p) , (3.13) 


with = uj. Similarly, one transforms the RHS of Eqs. (USD and (jS3D using 
the above prescription 1) ~ 3). 


3.2 LO equations 

The LO equations are obtained by neglecting O {dxt^dp^ and higher 
contributions in the gradient expansion. To this order the transport 
equations then read: 

2p^dx.F{X-p) = tp{X-p)F{X-p)-tF{X-p)~Q{X-p) , (3.14) 
2p^^dx^Q{X-p) = 0. (3.15) 

We note that the compact form of the gradient expanded equation to LO 
is very similar to the exact equation for the statistical function (ESD- The 
main technical difference is that there are no integrals over the time history. 
Furthermore, to this order in the expansion the evolution equation for the 
spectral function p(X;p) becomes trivial. Accordingly, this approximation 
describes changes in the occupation number while neglecting its impact on 
the spectrum and vice versa. The LO equations provide the basis of almost 
all practical applications of transport and Boltzmann or kinetic equations. 

"'^The expression is actually correct to NNLO since the first correction is O ( (s^dxi^)^) ■ 



3.3 NLO equations 

To NLO, i.e. neglecting all O (^{dxi^dp^Y^ and higher contributions, the 
transport equations become substantially more involved: 

[2pf^dx. + {dx.M^ (X)) dp^] F {X-p) 

= t,{X-p)F{X-p)-tj.{X-p)~Q{X-p) 

+ i {tp {X;p) , ^ {S+ (X;p), F {X;p)] ,(3.16) 

[2p>^dx. + {dx.M^ (X)) dp^] ~q{X-p) 

= i{s,(X;p),p+(X;p)}^^ + i{s+(X;p),p(X;p)}^^ . (3.17) 
Here we have introduced the Poisson brackets 
{/(^;p) ;^(^;p)}^^ 

= (dpJ{X;p)^ {dx^g{X]p)) - (dxJ{X;p)^ {dpJ{X;p)) .(3.18) 

The somewhat tedious NLO expressions are more or less not applied in 
practice. The rapidly increasing complexity of transport equations at even 
higher order seems forbiddingly difficult. In this case one may better 
employ the comparably compact full expressions without gradient expansion 
anyways. 

4 Comparison of LO, NLO and “full” results 

Nonequilibrium dynamics requires the specification of an initial state. While 
the corresponding initial conditions for the spectral function are governed 
by the commutation relations (E3D, the statistical function F{x,y) and first 
derivatives at = |/° = 0 have to be specified. Following Ref. m , we 
consider a spatially homogeneous situation with an initially high occupation 
number of modes moving in a narrow momentum range around p^ = p\\ = 
±Pts along the “beam direction”. The occupation numbers for modes with 
momenta perpendicular to this direction, p\+p\ = p\., are small or vanishing. 
The situation is reminiscent of some aspects of the anisotropic initial stage in 
the central region of two colliding wave packets.® Of course, a peaked initial 

® Other interesting scenarios include “color-glass”-type initial conditions with 
distributions exp(—peaked around ps = 0 with “saturation” momentum Qs- 
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particle number distribution is not very specific and is thought to exhibit 
characteristic properties of nonequilibrium dynamics for a large variety of 
physical situations. 

The initial distributions do not correspond to thermal equilibrium, which 
is isotropic in momenta. Interactions will lead to the isotropization of the 
distribution, which finally approaches thermal equilibrium. An example of 
the earlier stages of such an evolution is shown in Fig. |21 at various times 
in units of the renormalized thermal mass rriR. Shown is the distribution as 
a function of p± (vertical) and p\\ (horizontal), where bright (dark) regions 
correspond to high (low) occupation numbers. One observes that for times 
exceeding some characteristic isotropization time of about ~ lOO/m/j the 
distribution starts to become rather independent of the momentum direction. 
For a time of 150/mj? the figure shows already an almost perfectly isotropic 
situation. It has been shown in Ref. ffH that the characteristic time scale 
for isotropization is well described by the relaxation or damping time 

^damp - , (4.1) 

where the imaginary part of the thermal equilibrium self-energy —is 
evaluated for on-shell frequency for momentum p^s- This time scale 
agrees with the characteristic time for the effective loss of details about 
the initial conditions jH IH]. The time tdamp for partial memory loss is an 
important characteristic time scale in nonequilibrium dynamics, which is 
very different from the thermal equilibration time. The distribution is still far 
from equilibrium and the approach to a thermal Bose-Einstein distribution 
takes substantially longer than tdamp- For the employed initial conditions 
and coupling strength the characteristic thermal equilibration times exceed 
t dam p by more than two orders of magnitude [TT] . 

More precisely, we investigate a class of spatially homogeneous initial 
conditions parametrized as 


F(a;°,|/°;^|^o 0=0 = 


no(^ -h 1/2 


iOri 


(4.2) 


with d^oF{x^,0-,p)\^o=Q = 0, d^odyoF{x°,y^-,p)\^o=yO=Q = [no{p) 1/2] Up for 
ojp = ^P± + Mq = M^(x° = 0). The initial distribution 

function UqIp) is peaked around the “tsunami” momentum p^g with amplitude 
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Figure 2: Snapshots at times rriR = 10, 50, 100 and 150 of the initially 
anisotropic occupation number distribution as a function of p±_ (vertical) and p\\ 
(horizontal). Bright (dark) regions correspond to high (low) occupation numbers. 
The shown resolution was achieved by using a 64^ spatial lattice. 

A and width cr: 

no(^ = A exp |-^ [pI + (h I - Pts)^] | • (4.3) 

We compute the nonequilibrium dynamics for a rather weak coupling of the 
^2<j)4_interaction. For Fig. |21 we use = 0.25 and an initial distribution 
peaked at pts = 0.8 m^i with A = 20 and a = 0.32 ru/j. For these parameters 
we obtain tdamp — 100/mij according to Eq. (I4.1|l . We will consider a range 
of different couplings and amplitudes below. 

The question is whether transport or kinetic equations can be used 
to quantitatively describe the early-time behavior before fdamp, which is 
necessary if their application to the problem of fast apparent thermalization 
explained in Sec. ^is viable. This is discussed in the following comparing to 
LO or NLO transport equations. 

4.1 Isotropization 

A characteristic anisotropy measure can be chosen as CH 

AF(X°; uj, Pts) = F{X°; cn, p) |p^=o,p||=pts “ P) |px=pt.,P||=o , (4.4) 

which vanishes for the case of an isotropic correlator since the latter depends 
on = p'j_ + p|. Fig. 01 shows the time evolution of the on-shell derivative 

dxoAF{X°;u{pts):Pts) HH- The thick curve represents the full result, which 
is obtained from solving the evolution equations (Ell) and (EH) for the self¬ 
energies EH and (ITTUl) without any gradient expansion. For comparison, 
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Xo [mR-'] 

Figure 3: The solid line shows the on-shell dxo^F in Wigner coordinates, as 
computed from the solution of Eqs. without a gradient expansion. 

The dashed line represents the same quantity using the expression from the LO 
gradient expansion according to Eq. (IXTH) . If the gradient expansion to lowest 
order is correct, then both lines have to agree. Eor this quantity one observes 
that the lowest-order gradient expansion becomes valid on a time scale of about 
the isotropization or damping time tdamp — lOO/m-ij. The inset shows the same 
results on a logarithmic scale to make the small on-shell dxoAg visible. Its full 
result is given by the solid lower line, whereas the LO contribution is zero. 


we evaluate the same quantity using the LO gradient expansion according 
to Eq. ()3.14ll . For this we evaluate the RHS of Eq. ()3.14D using the full 
result for F{X^-,uj,p) and g{X^-,uj,p). If the gradient expansion to lowest 
order is correct, then both results have to agree. Indeed, one observes from 
Fig. I21 employing the same parameters as for Fig. |2l that the curves for 
the LO (dashed) and the full result agree well at sufficiently late times. 
However, they only agree after some time, which is rather well determined 
by the characteristic time scale ~ Oamp for the effective loss of details 
about the initial conditions. The latter is given by the thermal equilibrium 
estimate iH). Similarly, one can observe this time scale from the decay 
of the unequal-time spectral function in thermal equilibrium, as 

explained below. The latter measures the characteristic decay of correlations 
at time t with the initial state. The decay time for coincides rather well 
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Figure 4: Same as in Fig. |31 but for amplitudes A = 2,3,4, i.e. smaller 
occupation numbers according to Eq. dSl), and coupling = 0.5. The lower 
solid/dashed/dotted curves represent the full result, while the upper (thick) ones 
show the respective LO approximation. 

with the time for the onset of validity of the LO result. 

Using smaller amplitudes A (with a and pts kept hxed) one observes a 
similar picture. In Fig. H] we compare runs for A = 2, 3 and 4 with = 0.5. 
For instance, for A = ?> the LO order and the full result approach each 
other rather closely after a time of about 500/m/j. This time is rather well 
described by the characteristic decay time of the unequal-time correlator 
as shown in the inset of Fig. El Of course, the precise notion of a time after 
which a transport description holds depends on the definition and prescribed 
accuracy. We hnd for A = 2,3,4 that the damping time according to Eq. (ED 
is Oamp = 270,183,141m^^. For A = 3 this time may be compared to the 
inset of Fig. El which shows that only after, say, a time of about the 

correlator is small enough that initial-time effects play no important role. 

One may expect that the failure of the LO transport equations to 
describe the dynamics before isotropization completes is due to substantial 
NLO contributions. However, this turns out not to be the case. Fig. El 
compares LO, NLO and full result for the time evolution of the on-shell 
dxoAF{X^-,uj{pts),Pts) for H = 3 and = 0.5. The solid line represents 
the result as computed from the solution of Eq. (ESD without a gradient 
expansion. As before, the dashed line shows the same quantity according 
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Figure 5: Comparison of LO, NLO and full result for the on-shell dxoAF with 
A = 3 and g'^ = 0.5. For this quantity LO and NLO results do not differ much, 
while there is a substantial deviation from the full result at early times. The 
inset displays the real-time equilibrium spectral function for the considered mode, 
which exhibits the characteristic damping time for (partial) memory loss. The 
latter agrees rather well with the time for the onset of validity of a transport 
description. 


to the LO expression given by the RHS of Eq. (ITTl) using the full results 
F{X^;uj,p} and g{X^;uj,p}. Similarly, we evaluate the RHS of the NLO 
expression given by the RHS of Eq. dnni).® The result corresponds to the 
dotted line in Fig. El which turns out to be rather close to the LO curve. 

For the quantity dxoAF{X^; ^{Pts),Pts) the NLO corrections remain small 
even for rather strong couplings. This is demonstrated in Fig. El for g'^ = 1 
and same initial conditions as for Fig. El The convergence of the gradient 
expansion itself is apparently quite good here even for couplings of order 
one. However, when studying the importance of NLO corrections one has to 
distinguish between quantities whose vanishing corresponds to an enhanced 
symmetry from more generic quantities. Once an ensemble is isotropic it 
always remains so at later times even though the dynamics may still be far 
from equilibrium. The vanishing of AF obviously corresponds to an enhanced 
symmetry, i.e. rotation symmetry. In the following we consider more generic 

®More precisely, we employ the Poisson-Brackets where the djf-derivatives are 

calculated from the LO expressions (imi and EII3, which is correct to NLO. 
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Figure 6: Same as in Fig. El but for = 1. Even for this stronger coupling the 
NLO corrections remain small, while there is a substantial deviation from the full 
result at early times. 

cases. 

4.2 Thermal equilibration 

It has been demonstrated in Ref. m that in a scalar theory isotropization 
happens much faster than the approach to thermal equilibrium. Only a 
subclass of those processes that lead to thermalization are actually required 
for isotropization: “ordinary” 2 2 scattering processes are sufficient in 

order to isotropize a system, while global particle number changing processes 
are crucial to approach thermal equilibrium. 

Above we have investigated a quantity measuring an anisotropy 
(SjsfoAF), and which, therefore, approaches zero rather quickly after 
tdamp- Consequently, it is not affected by the relevant late-time processes 
which finally lead to a Bose-Einstein distribution. The accuracy of the 
gradient expansion for the late-time behavior can be investigated considering 
dxoF{X^]u},p) instead, which does not vanish for an isotropic ensemble. 
Fig. 13 shows the time evolution for this quantity on shell at zero momentum 
for the same parameters as for Fig. El The solid lines give the full result for 
dxoF (upper curve) and for dxog (lower curve). This is compared to the LO 
(dashed) and NLO (dotted) results evaluated in the same way as explained 
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Figure 7: Comparison of LO, NLO and full result for the on-shell dxoF at zero 
momentum with A = 3 and = 0.5. The lower curves show the full and NLO 
result for dxog, to lowest order there is no contribution. The LO transport equation 
fails to describe the full results until rather late times. In contrast, taking into 
account the substantial NLO corrections the gradient expansion becomes quite 
accurate for times larger than about the characteristic damping time for partial 
memory loss. 

above. 

One observes from Fig. [7|that the LO and the full result approach each 
other only at a time much later than it happens for the anisotropy quantity 
AF discussed in Sec. I4.1I However, one also observes from Fig. [3 that here 
the NLO corrections are substantial. The NLO results are much closer to 
the full result for times larger than about ~ tdamp- We emphasize that the 
employed coupling = 0.5 is not very large. The situation is qualitatively 
not very different also for couplings of order one, which is exemplified in 
Fig.IHlfor the same parameters as for Fig.lHl One observes, however, that the 
deviations of the NLO from the full result are stronger even at rather late 
times in this case. 

We conclude that quantities which are sensitive to the late-time dynamics 
can receive important corrections at NLO in the gradient expansion. This 
finding is, of course, not surprising for strong couplings. However, sizeable 
NLO corrections can be obtained already for couplings of about ~ 1/4. 
Nevertheless, we observe that for not too early times the NLO result can 
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Figure 8: Same as in Fig. Obut for = 1. 


get rather close to the full result for < 1. The agreement of the (NLO) 
transport description and full results start at about the same time as for 
the isotropization case discussed above. Times shorter than about Oamp are 
clearly beyond the range of validity of transport equations even for weak 
couplings. This limitation is apparently not a consequence of a gradient 
expansion to low orders. The latter seems to converge rather well for weak 
enough couplings. 

The failure of transport equations to describe early-time dynamics seems 
mainly based on the fact that the underlying assumptions 1) and 3) of 
Sec. 13.II are not controlled by a small parameter. The validity of the standard 
assumptions 1) and 3) both require an efficient loss of correlations with 
the initial state. The characteristic time scale for this loss is given by the 
inverse damping rate of the unequal-time correlator or, similarly, 

by (<F(f,a;)<F(0,^)). The latter describes the correlations with the initial 
time to = 0. A vanishing (<h(f, a;)<F(0, ^)) corresponds to an effective loss of 
information about details of the initial conditions. The characteristic time 
scale for this loss is well described by fdamp as dehned in Eq. dn} nn.^ 

^One may ask whether initial correlations beyond the considered Gaussian initial 
conditions could change this situation. However, such initial correlations would be beyond 
a typical transport description based on two-point correlation functions or effective particle 
numbers. 
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5 Off-shell transport 

The gradient expansion underlying transport equations is an expansion in 
the number of derivatives with respect to the center coordinates ()3.5|) and 
derivatives with respect to momenta associated to the relative coordinates 
(uni). The LO equations are obtained by neglecting O {dxi^dp^) and higher 
contributions in the gradient expansion, the NLO equations by neglecting 
O {{dxi^dp^^'^ etc. as described in Sec. El This expansion is constructed 
to work well for on-shell quantities, which are evaluated at the peak of the 
spectral function p(X°;a;,p) as a function of uj for given p. Away from this 
peak the contributions coming from derivatives with respect to momenta are 
typically not small if the spectral function is sufficiently narrow. The latter 
is always the case for weak enough couplings. Therefore, one may expect 
transport equations to fail for off-shell quantities even for weak couplings. 

In order to quantify this expectation, we have in mind a nonequilibrium 
situation for a given time A° larger than ~ tdamp but short compared 
to characteristic thermal equilibration times. In this case the occupation 
number can still show substantial deviations from a thermal distribution 
though certain properties, such as a fluctuation-dissipation relation 

F{uJ,p) = (^{uj)+ ^ , (5.1) 

with a non-thermal distribution n(a;), can become approximately valid |2I|. 
We assume that the spectral function may be approximated by its equilibrium 
form 


g{uj,^ = - , -2, (5.2) 

(^_^2+^ + M2 + 1S+(o;,p)) +(iS,(o;,p)) 

while the non-thermal distribution will be parameterized as (cf. also 
Eq. dOD) 

n(a;) = Aexp |-^(a;-a;o)^| . (5.3) 

We employ Eq. m with Eqs. dSH and ()5.3|) as a convenient class 
of approximate parametrizations of typical nonequilibrium situations at 
intermediate times as have been studied in the previous sections. With the 
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(a)-m)/y (co-m)/y 

Figure 9: The different curves of each graph show the different contributions 
FSg — (solid line), dxoM’^d^^F (dashed) and Q+}pb + ^{^+, F}pb 

(dotted). The first corresponds to the LO contribution in the gradient expansion, 
while the latter two add up to the NLO contribution. 


self-energies dZSI) and (j 2.1 Of) in Fourier-space (cf. See. 13.1^1 these form a closed 
set of self-consistent equations for F{uj,p) and g{uj,p). 

From these solutions we evaluate the LO and NLO expressions for 
dxoF {X^;uj,p) according to Eqs. (I3.14|l and (Id.lbll in the very same way 
as described in Sec. 0] The results are shown in Fig. El for various couplings 
as a function of (a; — for p = 0 with A = 1 and a = mp. Here 7 

corresponds to the half width of the spectral function. The different curves 
of each graph show the different contributions [FUg — pSi?]/ 2 u; (solid line), 

dxoM'^dt^F/2uj {dotted) scad {'Zp, 0 +}pb + F}pb /4a; (dashed) to the 

transport equation dnn). The hrst corresponds to the LO contribution in 
the gradient expansion, while the latter two add up to the NLO contribution. 
One observes that the LO contribution has its maximum value on shell, i.e. 
for uj = mp. For small enough couplings we clearly End that the on-shell 
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NLO corrections are small compared to LO ones. 

In contrast, Fig. El shows that for small couplings the NLO corrections 
have their maximum contributions off shell. They easily exceed the LO terms 
away from uj ~ rriR and here a gradient expansion to low orders cannot be 
expected to converge well. When the coupling grows to about g'^ ~ 1/4 the 
situation changes qualitatively. Here the Poisson-brackets contribution at 
NLO becomes comparable to the LO terms on shell. Though one expects 
transport equations to become unreliable for sufficiently large couplings, one 
observes for our example that NLO starts to become comparable to LO 
already for rather weak interactions. We emphasize that the normalization 
^gg (define our coupling reflects rather well the separation between 
a weakly -C 1) and strongly coupled regime {g"^ ^1). To see this we 
consider the ratio between damping rate and oscillation frequency of the 
unequal-time correlator which we denote as 27i'y/uj. Here 27i'y/uj 1 

corresponds to very weakly damped oscillations, while 27r7/a; S> 1 is 
overdamped. For comparison we display the corresponding values for this 
ratio in Fig. El One observes that they reflect rather well the coupling 
strength for the range of couplings considered here. 

6 Conclusions 

For the considered scalar g^^'^ theory we find that the earliest time for 
the validity of transport equations is well characterized by the standard 
relaxation or damping time tdamp given by Eq. (14.111 . We have done a 
series of comparisons for various couplings and conclude that for times 
sufficiently large compared to tdamp the gradient expansion seems to converge 
well. There are sizeable NLO corrections already for couplings of about 
g"^ ~ 1/4. Nevertheless, one observes that the NLO result can get rather 
close to the “full” result for g^ ^ 1. Times shorter than about tdamp seem 
clearly beyond the range of validity of transport equations even for weak 
couplings. This should make them unsuitable to discuss aspects of apparent 
early thermalization and stresses the need to employ proper equations such 
as (USD and (EH) for initial-value problems. 
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